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ABSTRACT 


An  analytical  method  was  developed  to  predict  the 
hydrodynamic  forces  and  moments  developed  on  a  body 
of  revolution  with  and  without  appendages.  The 
appendages  included  a  bridge  fairwater,  sternplanes, 
and  rudders.  It  was  assumed  that  the  body  was  either 
translating  with  an  angle  of  attack  of  up  to  20  degrees 
or  rotating  with  an  angular  velocity.  Calculations 
were  performed  by  using  the  analytical  method  for  four 
different  bodies  of  revolutions  and  comparisons  were 
made  with  available  experimental  data.  Potential  flow 
was  used  to  determine  the  inviscid  contribution  to  the 
hydrodynamic  forces  and  moments.  The  viscous  contri¬ 
bution  to  the  hydrodynamic  forces  and  moments  was 
determined  by  assuming  that  there  was  no  flow  separa¬ 
tion  due  to  the  axial  flow,  and  that  the  flow  separa¬ 
tion  on  the  leeward  side  of  body  due  to  the  crossflow 
was  independent  of  the  offsets  of  the  body,  the  angle 
of  attack,  the  angular  velocity,  and  the  forward  speed. 
It  was  assumed  that  the  angle  at  which  the  flow  sepa¬ 
rates,  as  measured  from  the  stagnation  point,  was  110 
degrees  at  the  stern  and  170  degrees  at  the  bow  for 
the  body  at  angle  of  attack  and  160  degrees  for  the 
body  with  an  angular  velocity. 


ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Submarine  Technology  Program  (STP)  Office 
of  the  Defense  Advanced  Research  Projects  Agency  (DARPA)  under  Task  Area  SI 974- 
030,  Program  Element  63569N.  The  David  Taylor  Research  Center  Work  Unit  Number 
was  1540-002. 

INTRODUCTION 


The  accurate  prediction  of  the  hydrodynamic  forces  and  moments  developed 
on  the  hull  and  appendages  of  a  submerged  vehicle  is  required  for  determining 
its  stability,  control,  and  maneuvering  characteristics.  Various  analytical 
methods  have  been  developed  to  make  these  predictions,  including  those  of  com¬ 
putational  fluid  dynamics,  but  none  of  them  have  been  successful  at  the  present 
time.  This  is  because  it  is  particularly  difficult  to  accurately  determine  both 
the  distribution  of  the  velocity  over  the  surface  of  the  hull  and  the  location 
of  the  separation  lines.  A  limited  amount  of  experimental  data  exist  for 
various  bodies  of  revolution  at  zero  angle  of  attack  and  for  a  spheroid  at 
small  angles  of  attack. 

When  a  streamlined  body  of  revolution  is  moving  through  the  fluid  with  an 
angle  of  attack,  the  flow  separates  from  the  hull.  The  location  of  the  lines  of 
separation  has  an  important  effect  on  the  magnitudes  of  the  hydrodynamic  forces 
and  moments  developed  on  the  hull.  For  example,  once  there  is  flow  separation, 
the  form  (pressure)  drag  becomes  much  larger  than  the  frictional  drag. 
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Computational  fluid  dynamics  can  be  used  to  solve  the  Navier-Stokes  equa¬ 
tions  and,  hence,  calculate  the  velocity  distribution  for  the  flow  around  a 
body  of  revolution  at  small  angles  of  attack.  However,  at  the  present  time 
there  are  no  validated  methods  for  determining  the  location  of  the  separation 
lines.  Due  to  the  difficulty  of  using  computational  fluid  dynamics  for  calcu¬ 
lating  the  location  of  the  separation  lines,  an  approximate  method  has  been 
developed  to  compute  the  hydrodynamic  forces  and  moments  for  a  submerged 
vehicle. 

As  discussed  in  Reference  1,  the  method  assumes  that  the  total  hydrody¬ 
namic  force  and  moment  can  be  divided  into  inviscid  and  viscous  contributions. 
Inviscid  flow  is  assumed  for  determining  the  velocity  and  pressure  distribu¬ 
tions  over  a  particular  cross  section  of  the  hull  from  the  stagnation  point  on 
the  windward  side  to  the  point  where  the  flow  separates.  After  the  flow 
separates  on  the  leeward  side,  the  pressure  is  assumed  to  be  oonstant. 

The  location  of  the  two-dimensional  crossflow  lines  of  separation  is  deter¬ 
mined  empirically.  The  angular  location  is  measured  from  the  windward  stagna¬ 
tion  point.  For  straightline  motion,  the  angle  is  selected  to  be  170  degrees 
at  the  bow  and  110  degrees  at  the  stern. 

The  viscous  contribution  to  the  hydrodynamic  forces  and  moments  is  deter¬ 
mined  by  solving  the  two-dimensional  boundary  layer  equation  for  a  body  of 
revolution.  It  is  assumed  that  there  is  no  flow  separation  in  the  axial  direc¬ 
tion.  The  forces  on  the  control  surfaces  are  computed  using  lifting  line 
theory. 

This  report  describes  the  method  for  calculating  the  inviscid  and  viscous 
contributions  to  the  hydrodynamic  forces  and  moments  on  the  hull,  describes 
the  method  for  calculating  the  forces  on  the  control  surfaces,  and  discusses  the 
results  of  the  calculations. 


INVISCID  FORCE  COMPUTATION 

The  coordinate  system,  oxyz  moves  at  a  speed  of  U  ,  which  is  the  mean  speed 
of  the  body  along  the  positive  ox-axis  {see  Figure  1)  for  the  straightline 
motion.  The  positive  oz-axis  is  always  directed  vertically  upwards.  The 
origin,  o  or  0  is  located  at  the  center  of  gravity  of  the  body.  The  OX-axis  is 
parallel  to  the  body  axis  and  the  OZ-axis  is  perpendicular  to  the  body  axis. 

When  the  body  is  undergoing  a  rotational  motion,  a  positive  drift  angle  is 
defined  as  shown  in  Figure  1. 

The  total  velocity  potential  can  be  expressed  as 


•KXry.z)  =  «()Q(x,y,z)  +  (|)(x,y,z) 


(1) 


where  (fi  is  disturbance  velocity  potential  due  to  the  body.  The  gradient  of  <(i 
is  given  as 


V(})^  =  ^"^o'  straightline  motion 

=  (Rn-y,  nx,  O),  for  rotational  motion 


(2) 
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The  disturbance  velocity  potential  t|)  can  be  expressed  in  the  integral  form  as 


4)(x,y,z)  ‘  ^  f!  G(P/Q)<y(Q)dS(q) 

s 


(3) 


where  P(x,y,2)  is  the  field  point,  Q(^ft»yo*z  )  the  source  point,  S  the  body 
surface,  a  the  unknown  strength  of  sources  and  sinks  distributed  on  the  body 
surface,  and  G  is  the  Green  function  which  is  given  as 


Q{P,Q)  -  I 


(4) 


where 


The  unknown  strength  of  sources  and  sinks,  cr  is  determined  by  the  body  boundary 
condition  which  is 


{^//  Q(P,Q)<y{Q)dS{Q)}  «  -7(|.^.n  (6) 

The  forces  acting  on  the  body  is  expressed  as 


?  ■  -  //  pndS  (7) 

s 

where  n  is  the  unit  normal  vector  at  the  body  surface  and  p,  the  pressure 
around  the  body,  is 


p  «  §  (8) 

where  p  is  the  density  of  the  fluid.  The  detailed  procedure  of  the  numerical 
solution  of  Equation  (6)  is  given  in  Reference  2.  The  body  surface  is  discre¬ 
tized  with  many  surface  elements.  The  unknown  strengths  of  sources  and  sinks 
are  assumed  constant  on  each  surface  element.  Then,  Equation  (6)  can  be 
expressed  as 


1  N 
-1-  I  a. 
^"i.1  ’ 


a* 


3n  M 


(9) 


When  the  kernel  integral  is  numerically  evaluated,  the  solution  of  Equation  (9) 
is  same  as  that  of  system  of  linear  equations. 
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VISCOUS  FORCE  COMPUTATION 


For  the  computation  of  the  axial  force  which  is  parallel  to  the  OX-axis 
the  method  developed  by  Young  (Reference  3)  is  used.  In  this  method  it  is 
assumed  there  is  no  separation  in  the  flow. 

The  total  drag  coefficient  is  expressed  as 


C 


A 


B 


.Vs. 2 


(10) 


where  S  is  the  surface  area  of  the  body.  The  momentum  area  of  the  wake  far 
downstream  is  defined  as 


"o  '  2"  /  0^  (1  - 


(11) 


where  u  is  the  velocity  in  the  wake  parallel  to  the  body  axis  and  U.  the 
Incoming  velocity  parallel  to  the  body  axis.  The  velocity  U.  can  be  expressed 
as 


U 
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UqCOSo,  for  motion  in  the  vertical  plane 

UqCOS/J,  for  motion  in  the  horizontal  plane  (12) 

RJlcos/J,  for  rotational  motion 


In  Equation  (10)  k.  is  the  momentum  area  at  the  tail  and  the  velocity  at  the 

edge  of  the  boundary  layer  at  the  tail.  In  order  to  compute  k  ,  the  laminar 
boundary  layer  equation  of  a  body  of  revolution  is  solved  up  to  the  transition 
point  and  then,  the  turbulent  boundary  layer  is  computed  up  to  the  tail.  The 
laminar  boundary  layer  equation  is  expressed  with  an  assumed  velocity  equation 
in  the  boundary  layer  as 


dx 


tm . 

u 


u 


+  2^U"g(X) 


(13) 


where 
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2  = 

X  S  U'2  =  2(^) 

'dx' 

f(X)  =  (4  -  fffx  .  3II5X2  .  ^X3)/C 

"  (ifs^  - 

■  <516  *  55651/'= 

C  .  3^  _i.  S  2 

315  315  ■  9072 


(14) 


In  Equstion  (13),  x  is  the  distance  measured  along  the  surface  of  the  body  from 
the  stagnation  point,  r  is  the  radius  of  the  body,  £  is  the  boundary  layer 
thickness,  U  is  the  velocity  at  the  edge  of  the  boundary  layer,  and  v  is  the 
kinematic  viscosity  of  the  fluid.  The  equation  of  the  momentum  area  for  the 
turbulent  layer  is  expressed  as 

^  +  ^(H  +  2)k  n  ^2ffr  (15) 


where  H(«1.4)  is  the  ratio  of  displacement  thickness  to  momentum  area  and  t  is 
the  skin  friction  at  the  body  surface,  When  Equation  (16)  is  solved  up  to°the 
tail,  the  total  drag  can  be  computed  using  Equation  (10).  The  details  of  the 
derivation  of  Equations  (13)  through  (15)  are  given  in  References  3  and  4. 

The  method  to  compute  normal  force  and  its  moment  is  slightly  different 
from  that  to  compute  axial  force.  The  sectional  crossflow  drag  is  integrated 
along  the  axis  of  the  body.  The  two-dimensional  flow  is  given  by 


U  ■  ZUgSine  (16) 

where  U2  is  the  two-dimensional  free  stream  velocity  (see  Figure  2).  The  velo¬ 
city  Ug  can  be  expressed  as 

U2  =  U^sino,  for  motion  in  the  vertical  plane 

«  UgSinp,  for  motion  in  the  hori2ontal  plane  (17) 

52y,  for  rotational  motion 
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Furthermore,  it  is  assumed  that  the  flow  separates  and  the  separation  point  must 
be  known  to  perform  the  computation.  The  crossflow  drag  is  expressed  as 


D 


2D  ' 


Dp  -Dp 


(18) 


where  Dp  is  the  frictional  drag  and  Dp  is  the  f:>rm  drag. 

For  a  two-dimensional  section,  the  lamin'  ‘joundary  layer  equation  becomes 


dx 


^  +  z2u"g(X) 


(19) 


and  the  skin  friction  is  given  by  (see  Reference  3) 


T 


0 


DVlk  ■»  12)U 
66 


(20) 


The  turbulent  boundary  layer  equation  is  the  same  as  Equation  (1S)  with 
omission  of  2ffr  on  the  right-hand  side.  The  frictional  drag  can  be  computed  by 
the  integration  of  skin  friction  t  along  the  body  axis  after  Equations  (15), 
(19),  and  (20)  are  solved  numerically. 

The  pressure  in  the  separated  region  is  assumed  constant  and  the  pressure 
coefficient  is  expressed  as 

Cp  =  1  “  (^)^  when  e  < 

(21) 

-  1  -  (^)^  when  0  >  e^p 

where  U  is  the  velocity  at  the  point  of  separation.  Then  the  form  drag  is 
computed^as  follows: 

2  "" 

Dp  =  pu  /  Cprcose  de  (22) 


where  0  is  the  angle  at  the  point  of  separation  as  shown  in  Figure  2. 

The  velocity  U  in  Equations  (13)  and  (19)  is  the  tangential  velocity  in  the 
axial  and  crossflow  direction,  respectively.  Even  though  the  tangential  velo¬ 
city  in  the  crossflow  direction  is  simpiy  expressed  with  Equation  (17),  this 
velocity  is  computed  around  the  body  section  with  application  of  potential 
theory  as  described  in  the  previous  section. 
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CORRECTION  TO  AXIAL  VELOCITY 


The  velocity  distribution  on  the  body  surface  can  be  computed  by  differen 
tiating  Equation  (1)  with  respect  to  x,  y,  and  z  as 


0  =  (u,  V,  w)  =  +  (j))  (23) 

The  tangential  velocity  on  the  body  surface  in  the  axial  direction  is 

U  =  ^-t^  (24) 

where  t.  is  the  unit  tangential  vector  in  the  axial  direction  and  it  can  be 
expressed  as 


t,  =  tg  X  n 


(26) 


where  t,  is  the  unit  tangential  vector  in  the  crossflow  direction  and  given 
in  Equation  (27).  The  velocity  of  the  airship  Akron  expressed  with  Equation 
(24)  is  plotted  in  Figure  3.  The  computed  velocity  agrees  very  well  with  the 
measured  velocity  up  to  96%  of  the  body  length,  This  velocity  becomes  zero  at 
the  stern.  However,  the  measured  velocity  shown  in  Reference  5  does  not  become 
zero  at  the  stern.  In  order  to  make  the  computed  velocity  close  to  the  measured 
one,  the  velocities  near  the  stern  are  extrapolated  as  shown  in  Figure  3.  The 
extrapolated  velocity  distribution  agrees  well  with  the  measured  velocity.  The 
extrapolation  scheme  for  the  velocity  at  the  stern  was  adopted  by  Cebeci  and 
Smith  for  the  solution  of  three-dimensional  boundary  layer  equations  (Reference 
6). 


TANGENTIAL  VELOCITY  IN  THE  CROSSFLOW  DIRECTION 
The  tangential  velocity  in  the  crossflow  direction  can  be  expressed  as 


U  =  V*t2 


(26) 


where  is  the  unit  tangential  vector  in  the  crossflow  direction  and  it  can 
be  expressed  as 
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tg  =  (cosflsina,  sine,  cosGcoso)  or 


(27) 


=  (sinesin/S,  sinecos/S,  -cose)  with  0  =  tan"\z/Y) 


The  computed  tangential  velocity  in  the  crossflow  direction  for  the  straight- 
line  motion  is  plotted  in  Figure  4  and  it  is  compared  with  Equation  (16). 

The  agreement  is  very  good.  This  indicates  that  for  the  straight! ine  motion 
the  tangential  velocity  in  the  crossflow  direction  can  be  computed  with 
Equation  (16).  The  direction  of  the  flow  is  always  in  one  direction.  The 
computed  tangential  velocity  in  the  crossflow  direction  for  the  rotational 
motion  is  plotted  in  Figure  6.  Compared  with  the  results  of  Equation  (16),  the 
computed  velocity  does,  not  agree  well.  Furthermore,  there  is  change  in  flow 
direction  depending  on  the  location  of  the  section.  The  location  of  the  maximum 
velocity  is  always  90  degrees  from  the  horizontal  plane. 


FORCE  COMPUTATION  FOR  THE  CONTROL  SURFACES 

To  compute  forces  and  moments  of  the  control  surfaces,  the  method  which  has 
been  developed  for  the  airplane  wing  by  Sivells  and  Neely  (Reference  7)  are 
used  herein.  In  this  method  the  distributed  lift  along  the  wing  span  is 
computed  with  the  lifting  line  theory.  The  sectional  lift  is  computed  at  the 
effective  angle  of  attack  which  can  be  expressed  as 


a.  »  a.  -  o. 
0  0  1 


(28) 


Where  o  is  the  geometric  angle  of  attack  and  Cj  the  induced  angle  of  attack 
(see  Figure  6) . 

To  compute  the  induced  angle  of  attack,  the  half  span  of  the  control  plane 
is  divided  as  follows 


y 


cos6. 


6=9,  18, 


81,  90 


(29) 


When  e  is  0,  the  section  is  located  at  the  tip  and  when  6  is  90,  the  section 
is  located  at  the  midspan. 

The  sectional  lift  is  approximated  as 


2C 


CT(y) 


A+1 .8' 


|0.6  . 


(30) 


where  A  (*b/C|^)  is  the  aspect  ratio,  c  is  the  mean  chord,  c^  is  the  chord  at 
y,  and  C,(a  )"^is  the  two-dimensional  lift  coefficient  for  a  . 

The  induced°ang1e  of  attack  is  expressed  as 


8 


(31) 


I 


d[CT(y)] 

Vi  -  y 


Once  the  induced  angle  of  attack  is  known,  the  induced  drag  is  computed  as 


®  C^tano^ 


(32) 


At  the  tip  of  the  control  surface,  the  induced  angle  of  attack  is  set  equal  to 
the  geometric  angle  of  attack  if  the  induced  angle  of  attack  is  calculated, 
using  Equation  (31),  to  be  larger  than  the  geometric  angle  of  attack.  Since 
the  calculated  sectional  lift  is  smaller  than  the  experimental  lift,  a  cor¬ 
relation  faction  of  0.8  is  applied  to  Equation  (31). 


LIFT  COMPUTATION  FOR  THE  TWO-DIMENSIONAL  WING  SECTIONS 

The  forces  on  a  two-dimensional  wing  section  can  be  computed  with  the 
vortex  panel  method  described  in  Reference  6.  The  sectional  contour  is  replaced 
with  a  two-dimensional  vortex  singularity  and  the  strength  of  this  singularity 
is  determined  by  the  body  boundary  condition  that  the  normal  velocity  at  the 
body  is  zero.  The  sectional  contour  is  discretized  with  m  straight  segments 
as  shown  in  Figure  7.  The  strength  of  vortex  singularity  varies  linearly  on 
each  segment.  The  velocity  potential  at  the  ith  control  point  (x,,y,)  can  be 
expressed  as 

y(s,)  .1  Vi-Vj 

tan  (33) 

where  >'(Sj)  is  the  strength  of  vortex  singularity  at  the  jth  segment  which  can 
be  expressed  as 


y(Sj) 


(y 


j+1 


-Vsj 


(34) 


S.  in  Equation  (34)  represents  the  length  of  the  jth  segment,  and  and  y.^. 
are  the  vortex  strengths  at  the  beginning  and  at  the  end  of  the  jth“'segmenti ' 
respectively.  The  (m-t'i)  values  of  y,  at  boundary  points  (X.,Y.)  are  unknown 
constants  to  be  determined  numerically.  The  boundary  condition  of  no  flow 
through  the  wing  requires  that  the  velocity  in  the  direction  of  the  outward  nor¬ 
mal  vector  be  vanishing  at  the  ith  control  point,  so  that 


an 


'l'(x^*y^)  *  0,  for  iel ,  2, 


(35) 
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The  unit  normal  vector  is  expressed  as 


n^  =  (-sine^,  cose^) 


(36) 


If  we  carry  out  the  involved  differentiation  and  integration  along  the  segments, 
we  obtain  the  followings 


where 


=  2ffUQSine^ 


for  i=1 ,  2,  ,  m 


(37) 


«  0.5FD  +  CQ  - 

C^2  "  °  +  FQ/(2Sj)  -  G(AC  +  DE)/Sj 
A  -  -[(x^  -  Xj)cos0j  +  (y^  -  Yj)sinej] 
B  *  (x^  -  Xj)^  +  (y^  -  Yj)^ 
c  ■  sin(e^  -  0j) 

D  *  cos(e^  -  ej) 

E  ■  (x^  -  Xj)sinej  -  (y^  -  Yj)cosej 
+  2AS, 

F  -  ln(1  +  -1— g - J) 


0 

p 

Q 


-1 

(x^  -  Xj)sin(e^ 


(X.  -  Xj)cos(e^ 


-  zej)  +  (y^ 

-  ZSj)  -  (y^ 


Yj)cos(e^  -  ZGj) 
Yj)sin(0^  -  aOj) 


(38) 
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When  -isj,  the  coefficients,  C,  E,  and  F  become  zeroj  and  C  .  and  C  -  in 
Equation(38)  have  simplified  values  as  follows 


Since  all  of  the  coefficients  in  Equation  (38)  are  dependent  on  the  shape  of  the 
sectional  contour,  the  coefficients  C  .  and  C  ,  can  be  evaluated  easily.  Once 
C  .  and  C  ,  are  numerically  evaluated?\he  unknown  strengths  of  vortex,  y.  can 
be^solved^using  Equation  (37).  However,  there  are  only  m  equations  in  Equation 
(37)  for  (m+l)  unknowns.  The  last  equation  which  supplements  Equation  (37)  can 
be  obtained  with  the  Kutta  condition  that  the  strength  of  vortex  at  the  trailing 
edge  be  zero  as  follows 


r,  *  r,*,  ■  0  (40) 

Once  the  vortex  strengths  are  numerically  solved,  the  velocity  components  can 
be  computed  through  differentiation  of  Equation  (33)  with  respect  to  x^  and  y.. 
The  pressure  on  the  contour  is  given  as 

P  -  f[Uo  -  |V«|/1^]  (41) 


The  force  and  moment  of  the  contour  can  be  computed  by  integration  of  the 
pressure  around  the  contour. 


INTERACTION  EFFECT  BETWEEN  BODY  AND  CONTROL  SURFACES 

As  indicated  in  Reference  9,  the  lift  developed  on  a  blunt-based  body  with 
fins  as  stern  appendages  can  be  calculated  by  using  slender  body  theory  to 
estimate  the  contributions  of  the  body  alone,  the  lift  on  the  fins  in  presence 
of  the  body,  and  the  lift  on  the  body  due  to  the  fins.  Alternatively,  the  lift 
on  the  fins  alone  can  be  estimated,  as  well  as  the  ratio  of  the  lift  on  the  fin 
body  combination  to  that  of  the  fins  alone.  In  order  to  estimate  the  lift  on 
the  fins  alone,  the  half-fins  are  joined  together. 

The  method  for  estimating  the  lift  developed  on  a  submarine-like  body  with 
fins  at  the  stern  is  discussed  in  Reference  13.  In  this  reference,  the  lift  on 
the  fins  alone  is  estimated  by  using  the  method  provided  in  Reference  16.  The 
contribution  of  the  lift  on  the  fins  in  presence  of  the  body  can  be  estimated 
by  multiplying  the  lift  developed  on  the  fins  alone  by  the  empirical  factor 
1  -  r/(2b),  where  r  is  the  maximum  radius  of  the  body  and  b  is  the  distance 
between  the  tip  chord  of  the  fin  and  the  centerline  of  the  body  (outreach)  for 
values  of  b/r  between  approximately  0.8  and  1.2.  Since  value  of  the  ratio  b/r 
is  approximately  1.0  for  the  analysis  performed  herein,  the  empirical  factor 
1  -  r/(2b)  would  be  approximately  0.5,  and  this  value  is  used  for  the  calcu¬ 
lations  herein. 
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DISCUSSION  OF  NUMERICAL  RESULTS 


The  results  of  calculations  that  have  been  made  using  the  analytical  method 
that  has  been  developed  are  discussed  herein.  The  empirical  method  for 
selecting  the  two-dimensional  separation  angles  in  the  crossflow  direction,  the 
geometric  characteristics  of  the  configurations  for  which  the  calculations  of 
the  hydrodynamic  forces  and  moments  were  performed,  and  the  results  of  the  cal¬ 
culations  for  the  unappended  hull,  the  control  surfaces,  and  the  partially 
appended  body  are  discussed  herein. 


TWO-DIMENSIONAL  SEPARATION  ANGLES  IN  CROSSFLOW  DIRECTION 

One  of  the  most  difficult  tasks  in  this  investigation  was  to  determine  the 
location  of  the  separation  line  on  the  surface  of  the  body.  The  hydrodynamic 
forces  and  moments  developed  on  the  body  can  not  be  predicted  without  knowing 
the  location  of  the  separation  line  accurately.  Even  though  the  velocity 
distribution  on  the  surface  of  the  body  is  computed  accurately  by  solving  the 
Navier-Stokes  equations,  there  is  poor  agreement  between  the  computations  and 
the  measured  values  of  the  hydrodynamic  forces  and  moments.  The  reason  for  the 
poor  agreement  in  the  force  and  moment  results  is  probably  the  inability  to 
predict  the  exact  location  of  the  separation  line  in  the  computation  process. 

While  the  velocity  distribution  on  the  surface  of  the  body  is  computed 
using  potential  flow  theory,  the  location  of  the  separation  line  must  be  accor¬ 
dingly  compatible  to  the  velocity  distribution.  The  location  of  the  separation 
line  is  a  function  of  the  Reynolds  number,  angle  of  attack  or  angle  of  drift, 
geometry  of  the  body,  and  the  turning  rate.  Figure  8  from  Reference  10  shows 
the  results  of  measurements  of  where  the  flow  separates  in  the  crossflow  plane 
of  a  spheroid  at  angle  of  attack  of  30  degrees.  The  figure  shows  that  the  angle 
at  which  separation  occurs  at  the  bow  is  larger  than  the  angle  at  the  stern. 

In  the  analysis  discussed  herein,  there  was  no  attempt  to  determine  the 
location  of  the  separation  line  analytically.  Instead,  the  hydrodynamic  forces 
and  moments  were  computed  using  an  arbitrary  location  for  the  separation  line, 
and  the  results  of  this  calculation  were  compared  with  experimental  data.  Based 
on  the  comparison,  the  location  for  the  separation  line  was  changed  and  another 
calculation  was  made.  After  one  or  more  additional  calculations,  the  computed 
results  for  forces  and  moments  and  the  experimental  results  agreed,  and  the 
location  of  the  separation  was  determined.  It  was  assumed  that  the  location  of 
the  separation  line  was  independent  of  the  Reynolds  number,  the  angle  of  attack, 
the  geometry  of  the  body,  and  the  turning  rate. 

The  fineness  ratios  of  the  nine  configurations  used  for  determining  the 
location  of  the  separation  for  straightline  motion  were  between  4.0  and  12.4. 

The  fineness  ratios  of  the  six  configurations  used  for  determining  the  location 
0*  the  separation  for  the  rotational  motions  were  between  4.0  and  11.0.  The 
values  for  the  separation  angles  that  were  selected  for  straightline  motion  were 
170  degrees  at  the  bow  and  110  degrees  at  the  stern.  The  values  selected  for 
rotational  motion  were  160  degrees  at  the  bow  and  110  degrees  at  the  stern. 

These  separation  angles  were  used  for  all  of  the  calculations  discussed  herein. 


12 


GEOMETRIC  CHARACTERISTICS  OF  THE  CONFIGURATIONS 

To  validate  the  present  method  the  forces  and  moments  of  four  configura¬ 
tions  were  computed  and  their  results  were  compared  with  test  results.  Their 
geometric  characteristics  are  given  in  Table  1.  For  a  spheroid  and  the  Akron, 
the  forces  and  moments  of  the  bare  hull  are  computed.  For  Model  4621  and 
the  DARPA  SUBOFF,  the  results  of  bare  hull  and  partially  appended  body  are  pre¬ 
sented.  The  total  number  of  discreti7ed  panels  for  the  body  surface  is2280  for 
all  of,tbe  models,  Forces  and  moments  are  non-dimensional ized  using  p\y\r /I 
and  pUqL®/2,  respectively.  ° 


RESULTS  OF  THE  CALCULATIONS  FOR  THE  UNAPPENDED  HULL 

To  validate  the  present  method  to  compute  the  axial  force  (Equation  (10)), 
the  computed  axial  forces  of  three  different  bodies  of  revolution  are  compared 
with  test  results  for  Reynolds  numbers  between  2  and  26  million.  The  results 
are  shown  in  Figure  6.  The  agreement  between  the  present  method  and  test 
(Reference  11)  is  very  good.  The  computed  potential  forces  should  be  zero. 
However,  small  values  of  forces  are  computed.  These  small  forces  are  caused  by 
the  numerical  inaccuracy.  The  body  surface  is  replaced  with  260  quadrilateral 
elements.  If  the  number  of  panels  is  increased,  the  potential  forces  will 
become  smaller. 

The  results  for  the  bare  hull  for  straightline  motion  are  given  in  Figures 
9  through  12.  Except  for  the  axial  forces,  the  computed  vertical  forces  and 
pitching  moments  of  spheroid  agree  very  well  with  test  results  (Reference  12) 
for  the  entire  range  of  angles  of  attack.  The  trend  of  the  measured  axial 
forces  is  different  from  those  for  other  configurations  as  will  be  seen  later. 

It  is  possible  that  there  might  be  some  errors  in  the  measurements.  The  poten¬ 
tial  axial  and  vertical  forces  become  zero  as  expected.  For  this  model  the 
numerical  inaccuracy  caused  by  discretization  of  the  forebody  cancels  out  that 
of  after  body. 

The  computed  results  of  vortical  forces  and  pitch  moments  of  the  Akron 
airship  agree  very  well  with  test  results  (Reference  5).  The  computed  axial 
forces  show  small  discrepancies  from  test  results.  The  discrepancies  are  caused 
by  the  fact  that  the  potential  axial  forces  are  not  completely  zero.  When  the 
angle  of  attack  is  zero,  is  computed  -0.0002.  If  the  number  of  panels  is 
increased  to  364,  C  is  computed  as  -0.0001.  The  error  analysis  is  given 
later  in  this  report. 

The  computed  results  for  the  forces  and  moments  for  Model  4621  and  the 
DARPA  SUBOFF  agree  very  well  with  test  results  as  shown  in  Figures  11  and  12. 

The  test  results  of  Model  4621  and  DARPA  SUBOFF  are  taken  from  References  13  and 
14,  respectively. 

The  results  for  the  bare  hull  for  the  rotational  motion  are  given  in 
Figures  13  and  14.  The  computed  results  of  lateral  forces  and  yaw  moments  of 
the  spheroid  agree  fairly  well  with  test  results.  The  measured  axial  forces  are 
not  available.  For  this  model  the  optimal  two-dimensional  separation  angles  are 
actually  found  as  160  and  120  degrees  when  the  drift  angle  is  nogative,  and  160 
and  100  degrees  when  the  drift  angle  is  positive. 

The  results  for  Model  4621  are  given  in  Figure  16  for  four  different 
turning  rates.  All  computed  forces  and  moments  agree  very  well  with  unpublished 
test  results.  However,  for  small  turning  rates  for  which  the  model  is  far  away 
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from  the  center  of  rotation,  the  computed  axial  forces  at  a  drift  angle  of  20 
degrees  are  larger  than  the  test  results,  while  for  turning  rates  larger  than 
0.39  the  computed  axial  forces  agree  very  well  with  test  results. 


RESULTS  OF  THE  CALCULATIONS  FOR  THE  CONTROL  SURFACES 

To  validate  the  present  method  of  computing  the  forces  and  moment  of  a 
finite  wing,  which  is  described  with  Equations  (28)  through  (32),  the  forces 
and  moments  of  a  finite  wing  with  three  different  aspect  ratios  are  computed  and 
the  results  are  given  in  Figure  15.  The  section  shape  of  the  wing  is  NACA  0015, 
its  taper  ratio  is  0.45,  and  the  sweep  angle  is  zero.  The  details  of  the 
geometry  and  test  results  ace  given  in  Reference  15.  Forces  and  moments  are 
non-dimensional i zed  with  pU^S  /2  and  pU^bS  /4,  respectively. 

The  computed  results  for^drag  and  lift  agree  well  with  test  results  for 
all  aspect  ratios.  The  computed  pitching  moments  with  respect  to  the  mean 
quarter-chord  agree  fairly  well  with  test  results  when  the  angle  of  attack  is 
less  than  12  degrees  and  when  the  aspect  ratios  are  1  and  2.  When  the  angle  of 
attack  is  larger  than  12  degrees,  the  computed  pitching  moments  do  not  agree 
well  with  the  teat  data.  When  the  aspect  ratio  is  3,  the  computed  pitching 
moments  agree  fairly  well  with  test  results  for  the  entire  range  of  angles  of 
attack.  An  improvement  in  the  theory  is  necessary  for  the  computation  of  the 
pitching  moment.  It  should  be  emphasized  that  a  small  error  in  predicting  the 
center  of  pressure  on  a  stern  control  surface  does  not  have  a  significant 
effect  on  the  estimate  of  the  pitching  moment  which  is  resolved  about  the 
center  of  buoyancy. 


RESULTS  OF  THE  CALCULATIONS  FOR  THE  PARTIALLY  APPENDED  BODY 

The  appendages  used  in  the  present  validation  are  listed  in  Figure  17. 
Appendages  3A,  11A,  3B,  3C  and  the  appendages  for  the  SUBOFF  are  four  identical 
cruciform  planes.  It  is  assumed  that  the  geometrical  angle  of  attack 
(Equation  (28))  for  the  sternplanes  or  rudders  for  the  straightline  motion  is 
the  same  as  the  angle  of  attack  or  drift  angle  of  the  body.  The  computed  forces 
and  moments  of  Model  4621  with  Appendages  3A,  11  A,  and  3B  agree  very  well  with 
the  test  results  for  the  entire  range  of  angles  of  attack  as  shown  in  Figure  18. 
In  the  computed  vertical  forces  and  pitching  moments,  there  is,  of  course,  no 
contribution  of  the  rudder  to  the  total  normal  forces  and  moments.  However,  in 
the  computed  axial  forces,  there  is  a  contribution  from  rudder. 

The  computed  results  for  the  SUBOFF  model  with  bridge  fairwater  and  stern 
appendage  egree  very  well  with  test  results  for  drift  angles  between  -20  and  20 
degrees  as  shown  in  Figures  19  and  20. 

The  computed  results  for  Model  4821  with  Appendage  3C  are  plotted  in  Figure 
21  for  three  different  turning  rates.  When  an  appended  body  is  undergoing  a 
rotational  motion,  the  geometric  angle  of  attack  of  the  rudder  in  Equation  (28) 
is  not  same  as  the  drift  angle.  For  Model  4621  it  has  been  found  that  the 
geometric  angles  of  attack  are  3.4,  9.5,  and  15.1  degrees  larger  than  the  angle 
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of  drift  at  the  center  of  rotation  for  turning  rates  0.142,  0.399,  and  0.639, 
respectively. 

As  shown  in  Figure  21,  the  computed  results  agree  well  with  test  results. 
Even  though  test  results  are  shown  for  drift  angles  are  between  -5  and  10 
degrees,  the  numerical  results  indicate  good  agreement  over  the  full  range 
of  higher  drift  angles. 


UNCERTAINTY  ANALYSIS  FOR  THE  NUMERICAL  PROCEDURES 

The  major  error  in  the  numerical  procedures  for  the  present  method  occurs 
in  the  numerical  evaluation  of  the  kernel  function  of  Equation  (9).  The  surface 
area  of  the  body  is  discretized,  the  source  and  sink  strengths  are  assumed 
constant  on  each  discretized  element,  and  the  double  integral  in  Equation  (9) 
is  evaluated  numerically  on  each  element.  For  the  present  presentation  of 
numerical  results,  the  number  of  discretized  elements  is  280  for  all  of  the  con¬ 
figurations.  The  larger  the  number  of  elements  is,  the  more  accurate  the 
numerical  results  are. 

To  estimate  the  errors  in  the  numerical  procedures,  the  computed  results 
are  compared  with  the  e»act  solution.  The  d'Alembert  paradox  says  that  the 
forces  around  a  body  in  potential  flow  become  zero  regardless  of  angle  of  attack 
(see  Reference  16).  Table  2  shows  the  potential  flow  results  for  the  axial 
and  normol  forces  of  the  Akron  airship  which  are  dependent  on  the  number  of  ele¬ 
ments.  As  shown  in  this  table,  the  computed  forces  becomes  smaller  and  smaller 
as  the  number  of  elements  increases.  When  the  angle  of  attack  is  zero,  the  nor¬ 
mal  forces  for  different  numbers  of  elements  are  all  zero,  While  the  flow  on 
the  body  surface  is  symmetric  with  respect  to  the  horizontal  plane  through  the 
axis  of  symmetry,  the  numerical  errors  for  the  upper  section  of  the  body  cancel 
out  the, errors  for  the  lower  section  of  the  body.  The  axial  forces  were 
calculated  to  be  non-zero,  and  these  forces  become  smaller  as  the  number  of  ele¬ 
ments  increases. 

When  the  angle  of  attack  is  21  degrees,  then  tho  normal  force  is,  of 
course,  non-zero.  In  this  case  the  flow  on  the  body  surface  is  not  symmetric 
with  respect  to  the  horizontal  plane  and  therefore,  the  normal  forces  ere  com¬ 
puted  as  non-zero  values. 

As  can  be  seen  in  Table  2,  as  the  number  of  elements  increases,  the  values 
calculated  for  and  C,  become  smaller  and  approaches  zero.  The  ratio  of  the 
values  of  and  to  the  values  of  Cj.  and  C-  calculated  for  280  elements  is 
tabulated.  ^As  can^be  seen  the  ratio  becomes  smaller. 


CONCLUSIONS 

The  analytical  method  to  predict  forces  and  moments  for  a  body  of  revolu¬ 
tion  with  and  without  appendages  has  been  developed  and  the  numerical  results 
for  several  configurations  show  very  good  agreements  with  test  results.  The 
present  method  can  be  applied  to  predict  the  in-plane  forces  and  moments  of  body 
of  revolution  with  and  without  appendages.  The  present  method  can  be  easily 
extended  to  the  computation  of  roll  moment  and  the  out-of-plane  forces  and 
moments. 
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RECOMMENDATIONS 


It  is  recommended  that  the  following  work  be  undertaken  to  upgrade  the 
method  developed  herein: 

1.  The  method  should  be  extended  to  include  the  prediction  of  the  forces  and 
moments  for  fully  appended  bodies,  and  the  results  should  be  applied  to  the 
analysis  of  the  stability  and  control  characteristics  of  submarines. 

2.  The  method  should  be  extended  to  include  free-surface  effects. 

3.  A  more  accurate  method  should  be  developed  for  determining  the  location  of 
the  lines  of  separation  in  order  to  improve  the  prediction  of  the  forces  and 
moments  at  large  angles  of  attack. 

4.  The  method  should  be  extended  to  predict  the  forces  and  moments  of  a  body 
of  revolution  undergoing  combined  pitching  and  yawing  angular  velocities. 
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SUBOf^F  with  stern  appendages 


Table  1  -  Geometric  characteristics  of  sample  models 


Model  name 

Length 

(ft) 

Beam 

(ft) 

Length/Beam 

Spheroid 

6.S6 

1.64 

4.00 

Akron 

19.62 

3.32 

5.91 

4621 

16.00 

2.04 

7.34 

SUBOFF 

14.29 

1.67 

8.67 

Table  2  -  Estimated  errors  in  the  numerical  procedures 
performed  for  the  Akron  airship 


0)  Angle  of  attack  ■  0  degrees 


Force 
coeff . 

Number  of 

elements 

?B0 

364 

420 

476 

O 

X 

-0.189x10'® 

-0.149X10“® 

-0.129X10"® 

-0.118X10“® 

V^'^X^280 

1.0 

0.80 

0.68 

0.62 

0.165x10'® 

0.194x10"® 

0.162x"® 

0.186x10"® 

V«'^2)280 

0 

0 

0 

0 

(2)  Angle  of 

attack  >  21  degrees 

-0.121x10'® 

-0.900x10“^ 

-0.746x10“^ 

-0.667x10"^ 

V^'^X^280 

1 .0 

0.74 

0.62 

0.54 

-0.297X10"® 

-0.223x10"® 

-0.187x10"® 

-0.164x10"® 

V^'^Z^280 

1.0 

0.75 

0.63 

0.55 
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